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Abstract. In tiny metallic nanostructures, quantum confinement and nonlocal response change the collective plasmonic 
behavior with important consequences for e.g. field-enhancement and extinction cross sections. We report on our most 
recent developments of a real- space formulation of an equation-of-motion that goes beyond the common local-response 
approximation and use of Ohm's law as the central constitutive equation. The electron gas is treated within a semi-classical 
hydrodynamic model with the emergence of a new intrinsic length scale. We briefly review the new governing wave equations 
and give examples of applying the nonlocal framework to calculation of extinction cross sections and field enhancement in 
isolated particles, dimers, and corrugated surfaces. 
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INTRODUCTION 

The interaction of light with the free electrons in noble 
metals has led to a range of novel plasmonic phenom- 
ena and a platform for a variety of new applications. 
In particular, nanofabrication technologies and chemi- 
cal synthesis are now allowing the plasmonics commu- 
nity to explore and manipulate light-matter interactions 
at sub- wavelength length scales, taking advantage of the 
spatially rapid oscillations of surface-plasmon polaritons 
and their ability to localize in very small metallic vol- 
umes and structures [1]. 

The understanding of the optical response of plas- 
monic structures has been successfully developed within 
the common framework of the local-response approxi- 
mation (LRA) with Ohm's law J(r) = cr(r)E(r) as the 
constitutive equation [2, 3]. However, the ability to fab- 
ricate and experimentally explore yet smaller metallic 
nanostructures has recently stimulated new theoretical 
developments aiming at an understanding of the nonlo- 
cal response and spatial dispersion [4, 5, 6, 7, 8, 9, 10, 
11, 12, 13, 14] as well as the onset of quantum con- 
finement effects and quantum tunneling as critical di- 
mensions are entering a new size regime where quan- 
tum mechanics is usually expected to become impor- 
tant [15, 16, 17, 18, 19]. While the commonly employed 
LRA is inherently a description without any intrinsic 
length scales [2, 3], the new developments naturally in- 
troduce new fundamental length scales associated with 



the quantum many-body wave dynamics of the electron 
gas. As a consequence, plasmon polaritons can not sus- 




FIGURE 1. Examples of metallic nanostructures, indicating 
characteristic length scales a for (a) isolated metal particles, (b) 
metal-particle dimers with gaps, (c) corrugated metal surfaces, 
(d) sharp metal tips, and (e) metal-nanowire metamaterials. 
While the local-response approximation is typically adequate 
for ka > 1, the nonlocal correction becomes important for 
ka<^ 1. 



tain spatial oscillations beyond the Fermi wave vector k F . 
For plasmon polaritons supported by an atomically thin 
sheet of graphene this becomes particularly clear [20] 
and there are similar cutoff effects in the homogenization 
of hyperbolic metamaterials [21]. However, even before 
reaching this cutoff we anticipate important nonlocal cor- 
rections in noble-metal nanostructures with critical di- 
mension approaching the 1-5 nanometer regime. 

For light interaction with arbitrarily shaped plasmonic 
structures, we review the first consistent real- space at- 
tempt of going beyond the LRA by a semi-classical hy- 
drodynamic equation-of-motion for the jelly response of 
the electron gas [8], while restraining ourselves from 
the sub-nanometer regime where the electron gas is an- 
ticipated to exhibit its full quantum wave nature, in- 
cluding quantum tunneling and quantum confinement 
effects. Having described the new nonlocal governing 
wave equations we proceed with some of our most re- 
cent analytical [8] and numerical explorations [10] of this 
framework in the context of field extinction and field en- 
hancement as well as implications for surface-enhanced 
Raman spectroscopy [11]. 

NONLOCAL RESPONSE THEORY 

We consider the interaction of light with metallic nanos- 
tructures in the linear regime, but with a nonlocal spatial 
response, i.e. 

V x V x E(r) = (f f J dr'e(ry)E(r'). (1) 

In the following we will focus on the electron plasma 
itself, for simplicity leaving out any interband effects. 
For later comparison, we note that within the LRA 
the two-point dielectric function simplifies to £(r,r') « 
£(r)5(r — r') with the usual Drude dielectric function 
e(r) = l + ia(r)/(£o(D) = 1 - G) 2 /[g)(g) + i/i)} where 
(O p is the plasma frequency, a is the Ohmic conductivity, 
and 1/t is the damping rate. 

HYDRODYNAMIC MODEL 

While nonlocal response, or spatial dispersion, is a con- 
sequence of the quantum many-body properties of the 
electron gas, we here limit ourselves to a semi-classical 
treatment [22, 23, 24, 25]. The usual equation-of-motion 
for an electron in an electrical field is extended to a hy- 
drodynamic one, including a pressure term originating 
from the quantum kinetics of the electron gas. By lin- 
earization, the plasmonic response is governed by the 
following set of coupled real- space differential equa- 



tions [8] 

VxVxE(r)=(f) 2 E(r) + /a)MoJ(r), (2a) 

^i)V[V. J(r)] + J(r) = a(r)E(r). (2b) 

Here, the V [V • J] = V x V x J + V 2 J correction to Ohm's 
law has a strength j3 2 = (3/5)v^ within the Thomas- 
Fermi model with vp being the Fermi velociy. Neglect- 
ing the quantum leakage of electrons (hard- wall confine- 
ment associated with a high workfunction) the boundary 
conditions for the current at the metal surface are particu- 
lar simple: the tangential component is unrestricted while 
the normal component vanishes as a consequence of the 
current continuity equation and the vanishing of all elec- 
tron wave functions right at the surface [8]. These cou- 
pled equations [Eq. (2)] form a natural starting point for 
numerical solutions of arbitrarily shaped metallic nanos- 
tructures [10, 26]. For numerical solutions we employ a 
finite-element method [10] and we have made our nu- 
merical implementation (an extension to the RF Module 
of COMSOL 4.2a) freely available [27]. Alternatively, 
for analytical progress one may proceed in the context of 
Eq. (1) and by eliminating the current from Eq. (2) we 
get 

£(r/) = S(r-0+G(r/)^ (3) 

with the dyadic Green's function defined by 

{^S 7 i) V [ V -] + 1 }G(r,r') = 5(r-r' ) . (4) 

Clearly, for j3 — >> the dyadic Green's function ap- 
proaches a Dirac delta function and we recover the usual 
LRA discussed below Eq. (1). Usually, this limit is con- 
sidered trivially fulfilled because vf <C c. However, do- 
ing a dimensional analysis the correction to the local- 
response Ohmic current has a strength (vf /c) 2 (ka)~ 2 
where k = co/c is the free-space wavenumber while a is 
a characteristic length scale, such as a nanoparticle diam- 
eter, a surface radius-of-curvature, or a gap distance sep- 
arating two nearby metallic structures, see Fig. 1. Thus, 
even though the electron dynamics is much slower than 
the speed of light the nonlocal response may nevertheless 
be an important correction if at the same time 

EXTINCTION CROSS SECTION 

The resonance in sub- wavelength (ka <C 1) isolated par- 
ticles [Panel (a) of Fig. 1] is commonly understood in 
the framework of polarization and the Clausius-Mossotti 
factor [2]. For spherical metallic particles this leads to 
a resonance in the extinction cross section at (Q p /V3, 
i.e. below the plasma frequency. Perhaps counterintu- 
itive, within the LRA the dipole resonance frequency is 



independent of the particle size. In the nonlocal treat- 
ment of the same problem, the resonance exhibits a size- 
dependent blue- shift, while there are no new resonances 
appearing below the plasma frequency [8]. On the other 
hand, above the plasma frequency a series of new reso- 
nances appear due to the excitation of confined longitudi- 
nal bulk plasmons. Already in the 1970' ties Lindau and 
Nilsson reported the excitation of bulk-plasmon reso- 
nances in ultra-thin Ag films [28, 29]. Here, we illustrate 
the same phenomena in Fig. 2 for Au spheres of radius 
a = 2 nm, where the scattering problem has been solved 
analytically, employing both quasistatic simplifications 
(V x E « 0) as well as by considering the fully retarded 
problem posed by Eq. (2). For the modeling of Au we 
use the following free-electron parameters: hcop = 8.812 
eV, H/t = 0.0752 eV, and v F = 1.39 x 10 6 m/s. Simi- 
lar blue-shifts where recently observed experimentally 
in EELS studies of Ag nanoparticles [30], though the 
results were interpreted in the context of a quantum- 
confinement model. Earlier elastic light- scattering exper- 
iments have indicated similar blue-shift effects [31]. 

The blue- shift is a natural consequence of nonlocal re- 
sponse; it requires a tiny extra fraction of excitation en- 
ergy to drive the jelly gas of electrons as compared to 
the LRA where each electron undergoes oscillations in- 
dependently of the other electrons in the gas, i.e. the in- 
ternal kinetic energy of the electron gas is not accounted 
for. As the radius is increased the nonlocal correction is 
reduced and the response of particles with a > lOnm is 
well accounted for by the LRA [8]. On the other hand, 
dimers of large particles with small gaps [Panel (b) of 
Fig. 1] can still support nonlocal effects even though the 
particles are far too big to exhibit any pronounced blue- 
shift themselves [4, 10] . In the LRA the induced charge is 
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FIGURE 3. Extinction cross section (7 ext as a function of 
frequency for TM-polarized light incident on a Au nanowire 
dimer with radius of 25 nm and a gap of 1 nm. The top panels 
show the local electrical field at the fundamental resonance for 
the LRA and the nonlocal model, respectively. 



strictly a surface charge. On the other hand, for nonlocal 
response the induced charge is smeared out over a finite 
volume near the surface with a characteristic length scale 
associated with the electron dynamics at the Fermi level. 
Obviously, this slight redistribution of charge influences 
both extinction and field enhancement in a quantitative 
way [10] and for kissing dimers the intrinsic length scale 
of the electron gas serves to resolve the otherwise diverg- 
ing behavior appearing within the LRA [12]. Fig. 3 illus- 
trates the extinction cross section for a dimer of two par- 
allel wires of radius 25 nm. Treated individually, their re- 
sponse is well accounted for by the LRA. However, when 
brought in close proximity to each other and separated 
only by a gap of 1 nm the nonlocal response quantita- 
tively changes the spectrum. In general, all spikes below 
the plasma frequency are blue- shifted and as a general 
trend, resonances at higher frequencies are more blue- 
shifted compared to those appearing at lower frequen- 
cies, as also observed in e.g. hollow nanowires [14]. 



FIELD ENHANCEMENT 



FIGURE 2. Extinction cross section c ex t as a function of 
frequency for TM-polarized light incident on a Au sphere in 
vacuum, comparing the fully retarded approach with the qua- 
sistatic approximation. 



Likewise to the extinction of dimers, the field enhance- 
ment is also influenced and in general the nonlocal re- 
sponse reduces the giant field enhancement supported by 
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FIGURE 4. Surface-averaged SERS enhancement factor for 
the case of Ag cylinder with radius 75 nm and with a 0. 1 nm 
radius-of-curvature smoothening of the contact point. The left 
inset illustrates a typical electric-field intensity of the funda- 
mental dipole mode, while the right inset shows the corre- 
sponding induced charge density. 



the tiny gap of the dimers. As in a classical capacitor, 
within the LRA the field enhancement in the gap diverges 
as the gap shrinks to zero [32]. On the other hand, non- 
local response prevents the field from diverging even in 
the limit of kissing particles [12]. This effect we have 
also studied for both dimers of cylindrical nanowires as 
well as for bow-tie type dimers [10]. Here, we turn to the 
problem of corrugated Ag surfaces [Panel (c) of Fig. 1] 
which also support huge field enhancement due to the 
deep grooves [33]. In our numerical examples we use 
Ag parameters from Ref. [34], treating interband effects 
within a single oscillator Lorentz model. In addition, 
we use vf = 1.3925 x 10 6 m/s appropriate for the free- 
electron response of silver. Within the LRA the SERS 

enhancement factor y(r) = |E(r) | 4 / |Eo | 4 would eventu- 
ally diverge when approaching the bottom of an infinitely 
well-defined groove [35]. On the contrary, due to the 
nonlocal response the field enhancement remains finite 
even in this extreme case [11]. In fact, even for rounded 
structures we find important quantitative differences be- 
tween the LRA and the nonlocal hydrodynamic model, 
as shown in Fig. 4 for the case of Ag half-cylinders 
resting shoulder-by- shoulder on an Ag substrate. In the 
particular example the cylinder radius is 75 nm and the 
contact point has been made smooth corresponding to a 
0.1 nm radius-of-curvature. Finally, Fig. 5 shows a typi- 
cal electrical-field intensity distribution for the excitation 
of an arbitrarily sharp Ag tip. Clearly, in the LRA such 




FIGURE 5. Typical electrical-field intensity for the excita- 
tion of an arbitrarily sharp Ag tip. The nonlocal response serves 
to smear out the diverging field of the otherwise singular geom- 
etry. 



a structure would exhibit a strong field singularity which 
is smeared out due to nonlocal response. 

CONCLUSION 

In conclusion, nonlocal response becomes noticeable in 
plasmonic structures with nanometer- sized features and 
critical dimensions. We have presented the first natural 
step in extending the common local-response approxi- 
mation to a semi-classical hydrodynamic treatment that 
takes the internal kinetics of the electron gas into ac- 
count. In contrast to the common treatment based on 
Ohm's law the new nonlocal framework offers an intrin- 
sic length scale inherent to the dynamics of the electrons 
at the Fermi level. As a consequence, field divergences 
are resolved even in geometries with abrupt and arbitrar- 
ily sharp changes in the surface topography. 
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